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We consider the electromagnetic field in the presence of polarizable point dipoles. In the corre¬ 
sponding effective Maxwell equation these dipoles are described by three dimensional delta func¬ 
tion potentials. We review the approaches handling these: the selfadjoint extension, regulariza¬ 
tion/renormalisation and the zero range potential methods. Their close interrelations are discussed 
in detail and compared with the electrostatic approach which drops the contributions from the self 
fields. 

For a homogeneous two dimensional lattice of dipoles we write down the complete solutions, which 
allow, for example, for an easy numerical treatment of the scattering of the electromagnetic field 
on the lattice or for investigating plasmons. Using these formulas, we consider the limiting case 
of vanishing lattice spacing, i.e., the transition to a continuous sheet. For a scalar field and for 
the TE polarization of the electromagnetic field this transition is smooth and results in the results 
known from the continuous sheet. Especially for the TE polarization, we reproduce the results 
known from the hydrodynamic model describing a two dimensional electron gas. For the TM 
polarization, for polarizability parallel and perpendicular to the lattice, in both cases, the transition 
is singular. For the parallel polarizability this is surprising and different from the hydrodynamic 
model. For perpendicular polarizability this is what was known in literature. We also investigate 
the case when the transition is done with dipoles described by smeared delta function, i.e., keeping 
a regularization. Here, for TM polarization for parallel polarizability, when subsequently doing the 
limit of vanishing lattice spacing, we reproduce the result known from the hydrodynamic model. In 
case of perpendicular polarizability we need an additional renormalization to reproduce the result 
obtained previously by stepping back from the dipole approximation. 

PACS numbers: 03.50.De 


I. INTRODUCTION 


Recently, in [1-3] the interaction of the electromagnetic field with a monoatomically this sheet polarizable perpen¬ 
dicularly was considered and different results were found. As an attempt to understand theses results we start from a 
two dimensional lattice of polarizable point dipoles in three dimensional space and reconsider the standard treatment. 
On this way, for a scalar field one arrives at the equation 


A-l-5^(5^(x-a„)j $(x) =0, (1) 

where the points On (n = (ni,n2,0), rii integer) form a plane two dimensional lattice of delta functions (’Dirac 
lattice’). 

Equations of such type, with delta function potential are known in different areas: 

1. The one-dimensional analog of eq. (1) with an equally spaced lattice, a„ = an, is the well known Kronig-Penney 
model (’Dirac comb’) [4], which now serves as a popular textbook example. 

2. In two and three dimensions, the delta functions in (1) are ill defined. For a single delta function this was 
observed in [5] and solved by the method of self adjoint extensions. 

3. In quantum mechanics, putting = 2mE/h^, eq. (1) is a Schrddinger equation with delta function potentials. 
It appeared as generalization of the Kronig-Penney model and provides an approximation for scattering of short 
wave lengths particles. It is known as method of zero-range potential (see [6] and references therein). Also, 
sometimes it is called point interaction. 

4. In quantum mechanics, in two or three dimensions, using methods known from quantum field theory, this 
equation is considered with some regularization supplemented by subsequent renormalization of the coupling 
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g. This was discussed in [7] and, more recently in [8], in application to delta functions on a manifold using 
a regularization in terms of separable potentials and application of renormaliation group techniques. A finite 
dimensional analog of this is known as Koster-Slater perturbation technique [9]. 

5. In classical electrodynamics with appropriate expression for g (see below), this equation describes the electric 
field in the presence of point dipoles, polarized by an applied field and by the field from all other dipoles. In 
this case, renormalization is equivalent to removing the self force. 

It must be mentioned that in each area the above mentioned methods are well developed, however, the communities 
using these, are quite disjunct. We find it useful to represent these methods together on one place demonstrating their 
close interrelation. We keep formulas as simple as possible to facilitate broader understanding and restrict ourselves 
to the examples relevant to the discussion of the polarizable sheets. 

The recent discussion of thin polarizable sheets started with [1, 10]. For a polarizability perpendicular to the sheet 
no interaction with the electromagnetic field was found. In [2] a response was found starting from a plane lattice of 
dipoles polarizable perpendicularly to the sheet. However, for a shrinking of the lattice spacing to zero, it became 
singular. In [3], also starting from a lattice, the limit of vanishing lattice spacing was taken in the equations. For 
polarizability parallel to the sheet, equations result which are equivalent to those of the hydrodynamic model, studied 
in [II]. These are of the type of eq.(l), however with a one dimensional delta function, (I(z), and have a well defined, 
unique solution. For perpendicular polarizability, however, an equation results, which has in addition a term with the 
second derivative of the delta function, S"{z). In order to match this term in the equation, the ansatz for the solution 
must itself contain a delta function. This is physically sound since the field $(x) describes in this case the normal 
component of the electric field of a double layer which is know to have a delta function on the sheet. However, in 
this case the product <I>(x)(I(z) is singular. As a way out, in [3] it was discussed that one needs to step back from the 
underlying dipole approximation. This way, a finite result was obtained, which is different from [I] and from [2]. It 
is an aim of the present paper, to discuss the limit of vanishing lattice spacing in a broader context and to gain more 
insight into the different limiting procedures involved. 

In electrodynamics, we consider the electric field in the presence of polarized point dipoles. The corresponding 
formulas are well known (see, e.g., [12]). Here we follow the notations used in [3]. Putting c = 1 and assuming 
harmonic time dependence ~ exp(—iwt), we have 2 equations, 

(—— A)E(x) = 47r(a;^ + V o V) ^ Pn^(x — a„), 

n 

Pn — Q;E(an)5 (2) 

where the first line is the equation for the electric field in the presence of the point dipoles Pn located at an, and the 
second line describes the polarization of one dipole by the electric field at the location of the dipole. In this way, the 
commonly used dipole approximation is assumed. In eq. (2), Vo V is the dyadic product of the gradients. In general, 
all vectors are denoted by bold letters, for instance 

X = (a;,y,2;) = (x||,z), X|| = (a;,y), (3) 

and similar for the momenta. The sheet is always in the (x,y)-plane and the lattice 

a„ = an, n = (ni,n2,0) (4) 

is quadratic with spacing a. In the second line in (2), d is the polarizability. It has the dimension of length^. We 
take it as a diagonal matrix and consider two cases, 

d = diag(a| |, a| |, 0) for in-plane polarizability, 

d = diag(0,0, as) for perpendicular polarizability. (5) 

If thinking of the dipoles as bound charges e with non relativistic motion, for a|| and as known expressions of the 
type 


a = 



( 6 ) 


with intrinsic frequency Wq hold in the simplest case. It is to be mentioned that the sign may change in dependence 
on the frequency. For wq = 0 we get the polarizability for free charges as assumed in the hydrodynamic model [II]. 
However, in this paper we do not use these details. 
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Next we insert the dipole moments from the second line in (2) into the first line and we obtain an equation like (1), 

— A — 47r(a;^ + V o V)q; ^ S(x — a„)^ E(x) = 0. (7) 

and get an effective Maxwell equationlt must be mentioned, that this equation is a direct consequence of the equations 
(2). However, it is common not to use this equation in electrodynamics, see for example Chapt. 17 in [12]. Below 
we will comment more on this point. Keeping aside for a moment the difficulties with the singularity, this is a vector 
equation. It can be diagonalized since the polarizations separate for the polarizabilities (5) and the geometry (4) into 
TE and TM modes. 


E(x) 


-dy\ 

dx Ete (x) 

0 J 


dxdz \ 

dydz i?TM(x) 

-All/ 


( 8 ) 


with A|| = d^+dy. For in-plane polarizability we get an equation for the TE mode where V o V does not contribute, 
whereas for the TM mode it becomes the in-plane Laplace operator, V o V —>■ A||. For perpendicular polarizability 
we have for parity reasons only the TM mode. In the following we denote this case be ’P’. We consider the equation 
for the normal component of the electric field, E 3 (x), and here we have to substitute V o V —>■ 9^. Together, in all 
three cases we get the modification 


- A + - an)J <l)(x) = 0. 

of equation (1) or (7) where for g and P one needs to substitute according to 


mode 

\ coupling 

field 

TE 

9 a||. 

P = -47ru;2, 

$(x) Ete(x), 

TM 

9 a||. 

P = -47r(u;^ -1- A||), 

4>(x) Etm(x), 

P 

g ^ as, 

P = -47r(u;2 -f92), 

$(x) Esix). 


(9) 


( 10 ) 


We mention that the TE mode differs from the scalar case only by the substitution of the coupling constant, while 
the other two cases have in addition derivatives in their couplings. 

The paper is organized as follows. In the next section we consider the equation with a single delta function and 
discuss in detail the singularity, the regularization and renormalization, which are used in section III for multiple 
centers. As an application we consider then the scattering on a two dimensional lattice. In section V we discuss the 
transition to a continuous sheet with different orders of the limits involved. 


II. APPROACHES FOR A SINGLE THREE DIMENSIONAL DELTA FUNCTION 

In this section we consider the approaches to eq. (1) for a single delta function located at the origin, a„ = 0, 

(-w2-A + 553(x))<i>(x)=0, (II) 

where d{x) = S(x)S{y)d{z) is a three dimensional delta function, x = {x,y,z) € and A is the three dimen¬ 
sional Laplace operator. As it stands, this equation is not well defined. In three (and in two) dimension this equation 
can be given a precise meaning by several the methods, which we review in this section. 


A. Self adjoint extension 

We begin with the method of self adjoint extensions. This was first done in [5], later generalized in [13] and quite 
a number of further places. The idea is, roughly speaking, to restrict the domain of the operator to regular solutions 
where it is only symmetric, and then to add a singular solution which makes the operator self adjoint. In doing so, 
a new parameter ose appears which can be called extension parameter and the original one, g, looses any meaning. 
In this way, any potential with support in one point is described. This solution appears only in the s-wave (orbital 



4 


momentum 1 = 0). Higher I are not affected (these ’do not feel’ the delta function) since any regular solutions behave 
~ for r —>■ 0. In this way, the solution has an expansion 

+ asE + 0{x)^ , (12) 

where c is an overall normalization constant and asE is the extension parameter. This extension is unique in the 
s-wave sector (see Sect. IIB). Another, equivalent form of writing, 



lim -asE + 
|xK0 


d|x| 


|x|$(x) = 0, or asE = lim In (|x|$(x)), 

|xKo d|x| 


is in terms of a boundary condition at the origin. 

The solution to eq. (11) for x ^ 0, i.e., outside the delta function, can be written in the form 


$(x) = + f 


piui\x\ 


(13) 


(14) 


It can be interpreted as incoming plane wave with wave vector k and an outgoing spherical wave with frequency oj 
(with uj = |k| from eq. (11)), which is centered around the delta function. In this setup, / is the scattering amplitude. 
Expanding this solution. 


$(x) = l + -^+za;/ + ... 

X 


(15) 


and comparing with eq. (12), we identify 


asE = T+*w, or, / =-(16) 

/ asE - 

which is a relation between the extension parameter asE and the scattering amplitude. In this way, the solution (14) 
is extended to the whole plane, i.e., including x = 0. 

The scattering amplitude / has the appropriate analytic properties. It is a meromorphic function with a single pole 
on the imaginary axis in w = i(— q;se), which for asE < 0 corresponds to a bound state. The normalized boundstate 
wave function is 


<(’bs(x) 



gasE|x| 

X 


(17) 


and (—ose) is the binding energy. 


B. The three dimensional delta function in the general theory of selfadjoint extensions 

In the last years Asorey, Munoz-Castaneda et al developed a new formalism to characterise the most general case 
of self adjoint extensions for operators relevant in mathematical physics in terms of an unitary operator U (see Refs. 
[14-18]). Here, we apply this formalism to a three dimensional delta function, which was not done so far. 

In general, if M is a d-dimensional manifold with boundary dM = S, the set of selfadjoint extensions of the 
Laplace operator on M is in one-to-one correspondence with the unitary group U (L^(M)). This correspondence 
allows to characterise each self adjoint extension by an unitary operator U. In particular, for a given unitary operator 
U GU (T^(M)), the corresponding self adjoint extension Afj is fully determined by the domain of functions V{Au) C 
L‘^{M) that satisfy the boundary condition 

/X V'Ie - * SnV’Is = C^(Ai V’Ie+* 5„V’ls)> (18) 

where the constant /x is a positive coupling constant with units of L~^ that characterises the interaction with the 
boundary S. When the size of the boundary is finite, the coupling constant /x can be set to unity without loss of 
generality since the spectrum of Ajj and its eigenfunctions only depend on the unitary operator U that characterises 
the self adjoint extension (see Refs. [15, 16]). In particular, the existence of bound states in the spectrum of A„ only 
depends on the matrix U and not in the value of /x. 
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In the case of a singularity concentrated on a point as it is the delta function, one must introduce a regularisation 
in order to define the singular point-like potential. For this reason it can not be set to one. Assuming that the 
singularity is placed at the origin, the operator U is defined on a sphere of radius e around the singularity. Note that 
in such a case the boundary S is a sphere of radius e and the physical space of the system is given by the condition 
r > e. Therefore we have i9„ = —dr and eq. (18) turns into 

Me i’\r=e + * dr'tp\r=e = Ue{^le “ * dri^\r=e)- (19) 

The radius e is here the regularising parameter that must be made 0 at the end. As it will be seen below, we need to 
allow the coupling constant fi to depend on the regularization parameter, fj, —>■ Therefore eq.(19) can be written 

as a boundary condition. 


lim (/Xe 'lp\r=e+idr'll^\r=e) = (20) 

When the singularity preserves spherical symmetry, the operator U, that characterises the corresponding self adjoint 
extension, is box-diagonal when the wave function is decomposed in spherical components. Therefore for each value of 
the orbital angular momentum L, we will have one boundary condition characterised by the a finite matrix of 

order 2L x 2L (see Refs. [19, 20]). Any quantum state tpki^) with defined energy can be decomposed in spherical 
coordinates as 


oo L 

ipk{^) = E E Rk,LM{r)yLM{^), (21) 

L=0 M=-L 

where are the spherical harmonics, and Rk,LM is the radial function that obeys the differential equation, 

— ('>'^-RLM{r)) + RkMir) = k‘^RLM{r). ( 22 ) 

T \ T ) 

The general solution to the radial differential equation is given by 

RLM{r) = ALM(k)jL(kr) + BLM{k)yL{kr), (23) 

where jL('r) and ykij) are the spherical Bessel functions. Notice that since the boundary condition in general does 
not require regularity of the radial function at the origin, the coefficient of yL^kr) can not be made equal to zero. 
Hence, for fixed L, the corresponding boundary conditions have the form 

L 

{e))YLMm (24) 

M=-L 


for any L. Since the spherical harmonics are orthonormal, these should be taken as the basis of a 2L vector space of 
boundary data for fixed L. The boundary data can be represented as two 2L-dimensional column vectors. 




( AieRfc.L-L(e) ± L_L(e) \ 


Me^feTM(e) ± iR'k.LM^^) 


\ yeRk,LL{e) ±iRk,LL(^) J 


(25) 


Now, the most general boundary condition that preserves the total angular momentum can be written as 

Ui G [/(2L + 1), L = 0,l,2,3,... . (26) 

This boundary condition is characterised by a set of matrices {UijUi G U{2L -|- l)}Lgz+. Notice, that when the 
boundary condition is defined by Ul = I 2 L-i-ix 2 L+i, one obtains the boundary condition R'k — 9 which is 

Neumann boundary condition. If Uz, = —I 2 L-i-ix 2 L+i gives rise to eRk^hhid) = 0 which must be satisfied for any 
e > 0 before taking the limit e —>■ 0 which happens if and only if Rk,LL{f) = 0 which is Dirichlet boundary condition. 
When taking the limit e —>■ 0 for the cases of Dirichlet and Neumann boundary conditions, the coefficient B^m in 
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equation (23) must be equal to zero which gives rise to functions regular at the origin. It is of note that when the 
point potential is spherically symmetric, the matrices must be diagonal and the boundary condition (26) becomes 

R'k,LMi^) = Me tan (0e(iy, lW)/2) i?fc^iM(e) (27) 


for each radial function Rk,LM- 

Within this formalism, it is also possible to define a multipole boundary condition by interpreting the L- 
decomposition of the boundary condition given above as a multipole decomposition of the point potential. Therefore, 
when Ul ^ —I 2 L+ix 2 L+i just for a given Lq > 0, then one should interpret the point potential as the short range 
approximation for the potential of a L-multipole (Lq = 0 would be the monopole or point charge, Lq = 1 the dipole, 
L = 2 the quadrupole, and so on). As it is well known the delta function potential represents the short range 
approximation of a point charge (this is the same as the zero range approximation discussed in the next subsection). 
Therefore, to represent the delta function potential as a boundary condition, only the monopole term in the multipole 
decomposition of the boundary condition should be non-trivial and = —I for all L > 0. 

To determine the L = 0 component we have to take into account that in this case Ul=o = e*® G U{1). Therefore 
the boundary condition for the L = 0 components reads 

Me7?fc,o(e) + iR'k,oi^) = e*^'Me7?fc.o(e) - iR'k, (28) 

To determine fie and in terms of e and the extension parameter asE introduced in (12), we need to write down 
condition (13) in a suitable manner. Previous to taking the limit |x| —>■ 0 in equation (14) we can write condition (13) 
over a small sphere of radius e, 

(1 - easE)Rk,oi>^) + e7?fc_o(*^) = 0 R'k^de) = —^ - Rk,o{<^)- (29) 

Provided ±1, we can write the condition (28) as 

775c.o(e) = Me tan {9J2) Rk,o{(^)- (30) 

Comparing the last two expressions we obtain the condition 

/ie tan (6>e/2) = , ( 31 ( 

Keeping in mind that the only dimensional parameter entering the boundary condition (18) is we obtain the 
relations 


= 1/e, tan {9J2) = easE - 1, 


(32) 


connecting the parameters of the self adjoint extension in subsection IIB and in the present subsection. 
It is of note that the renormalization group equations for the parameters /Xe and 9e are 


de 


-1/m?; 


de 


asE 

1 — tan^(0e) 


(33) 


This dynamical system also characterises the three dimensional delta function up to two integration constants. 

The generalisation of the definition of the delta function potential for an electromagnetic field can be achieved using 
the theory of self-adjoint extensions for the quadratic Yang-Mills operator around a point-like configuration (see Refs. 
[21-23]). 


C. Zero-range potential 

In quantum mechanics, eq. (11) with u> = \/2mE/h is a Schrodinger equation for a particle with mass m moving in 
a delta function potential. For a generic spherical symmetric potential V(r) in place of the delta function one considers 
the scattering of the particle off the potential and comes to a partial wave scattering amplitudes fi(k) (I = 0,1,2,, 
k = jkj). For small momenta k, or for large wave length of the scattered particle, these can be expanded, 

fi{k) = -aik'^‘ + ... 


(34) 
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(see, e.g., [24] chap. 11), where o; is the scattering length, which is a property of the potential. For small k, eq.(34) 
provides an approximation for the scattering amplitude. If restricting in addition to s-wave scattering, comparison 
with eq. (16) shows the relation 


ao = - — (35) 

asE 

between the extension parameter and the s-wave scattering length. For the s-wave, in the next term in expansion 

(34), 


fo{k) =-ao + ^rok'^ + ..., (36) 

the parameter tq is interpreted as range of the potential. Putting rp = 0 motivates the naming ’method of zero-range 
potential’. 

This method was actively used since the mid 1960ies, see, for example, [6] and also the more recent paper [25], 
which has an exhaustive bibliography. Some most recent applications are [26] and [27]. The equivalence to the method 
of self adjoint extension was shown in [28]. 


D. Regularization and renormalization 


In this subsection we consider another approach which follows a line of reasoning know in quantum field theory. 
Following [7], we start with a Fourier transform. 


*(x) = / ^ e- 0(p), 

and get from (11) 

(-w^ -k p^) (^(p) -k 54>(0) = 0 
with = k^. This equation can be easily solved. 


^(P) 


(27r)3(5(p-k) 


gm 


(37) 


(38) 


(39) 


assuming Sla; > 0. The first term is the homogeneous solution. Inserting (39) into (37) we get for x = 0 the equation 


where 


$(0) = l-g$(0)/(-w2-i0), 


I{z) 


d^p 1 

(27r)3 _|_ 2 


(40) 


(41) 


carries the divergence. This divergence is one way to see that the delta function in eq. (12) is ill defined in dimensions 
higher than one. At once, this integral gives the opportunity to introduce a regularization, i.e., to change the initial 
formal expression in a way that it becomes well defined. A convenient way is to change the dimension, 3 —>■ n, in 
the integration and, following the aim of dimensional regularization, to consider the analytic continuation of I{z) in 
n. Another way is to restrict the integration, |p| < A, by a momentum cut-off with A —>■ oo at the end. Thus one 
substitutes 


I{z)^l\z)= I 
|p|<A 


d^p 1 

(27r)3 + z 


A 

27r^ 




(42) 


and has explicitly the singularity for A —>■ oo. Rewriting eq. (40) with this regularization and solving for 5$(0), 


gm = ^ 

g 


ICO I 
477 


A 

2F 


1 


(43) 



one defines by 


1 _ 1 A 

9 t 9 27r2 


(44) 


a new, renormalized coupling constant gr which is supposed to be finite for A —>■ oo. This is on expense of the initial 
coupling constant g, which is now g = which goes to zero. Finally, inserting 54>(0) from (43) with 

(44) for A —>■ 00 into (39) we get the solution which after Fourier transform back becomes 


$(x) = 


ikx 


1 


i: + l|47r|x|- 


This is just the solution (16) with the identification 


j. —1 47r 

+ iw gr 

Dr 


(45) 


(46) 


In this way, one obtains after renormalization the same result as from self adjoint extension, eq. (16). The freedom 
in the renormalization (one can add to (44) any constant) corresponds to the freedom in the choice of the extension 
parameter asE- 

At this place it should be repeated that the initial parameter g in the process of renormalization looses its meaning 
completely. The new, renormalized parameter gr can be given a meaning by relating it via (46) with the scattering 
amplitude or, using (17) with a bound state level. 

The cut-off regularization introduced with eq. (42) is equivalent to the use of a regularized delta function (I(£,x) 
in eq. (11) which is non singular and has the property 


lim (j(£, x) = 5(x). (47) 

e—>0 


However, for any non zero e, one comes this way to an equation with a generic potential. Especially when considering 
several such potentials in an equation like eq. (1), variables do not separate and the equation becomes hard to 
investigate. There is a way to avoid this difficulty by considering the regularized equation in the form 


(-w^ - A)4>(a:)-hg(5i(e,x) J dy(j 2 (e, y)4’(y) = 0, 


(48) 


where we introduced two regularized delta functions which may be different from one another. 

Starting from here, we consider in parallel to the scalar case, eq. (1), also the electric field obeying eqs. (2), with 
s single dipole located at the origin, an = 0. We use this freedom in the choice of the regularized delta functions for 
accommodating the cases we are interested in. Thus we demand 

lim 5i(e, x) = PS{x.), lim 62 ( 6 , x) = 5(x). (49) 

e->0 e->0 

with P = 1 for the scalar case and with P either P —>■ —47r(w^ -I- V o V) (and g —>■ —47rd) or, according to formula 
(10), for the electromagnetic case. Clearly, for £ —>• 0, using (49), we get back to eq. (11) for the scalar case. For the 
electric field we substitute 4>(x) —>■ E(x) in (48) and get back eq. (7). 

Rewriting eq. (48) in the form 


J dyA:(x,y)$(y) = 0, 


(50) 


with the integral kernel 


K{x,y) = (-uj^ - A)d(x- y) + gSi(e,x)S 2 (s,y), 


it is seen that we have an equation with a separable potential. Defining the Green function by 


(51) 


dzK{x,z)G{z,y) =S{x-y), 


(52) 
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the explicit solution is 


(^(x, y) = G'o(x - y) - y dx'dy'Go{x - x.') 6 i{e, x')4 ^S 2 {s, y')Go{y' - y) 


with 


<Po = -+ [ dy:dy 62 {e,x)Go{x- y)Si{e,y) 

9 J 

and the free Green function obeying — A)Go(x — y) = (5(x — y). Its explicit form is 


Go(x- y) = 


dp 


^zpx 


3ZCl>|x| 


(27r)3 p2 — — zO 47r|x| 


(53) 


(54) 


(55) 


The solution (53) can be easily checked by inserting into (52). 

From (53), a solution corresponding to an incoming plane wave with wave vector k, can be obtained by 


$(x) = y dyd-z G(x, y)Go ^ (y - 


and with (53) it can be written in the form 

$(x) = - y dx'Go(x - x')(5i(e, x')(/)(C^ J dy' 52 {e, y') 

For £ —>■ 0 it turns into the solution (14) with the identification 




/ = lim 


-1 


e->o Aircfio 

This way, the equivalence of the separable regularization with the previously considered approaches is seen. 
An especially convenient choice for the regularized delta function is the heat kernel. 


(56) 


(57) 


(58) 


Ke{x) = 


exp - 


i-i) 


(47rs)^/^ ’ 

with its property lime_>o Kei^) = <5(x) and the pleasant formulas 

pOO 

Go(x) = / dse~‘'^^ Ks{x), 

Jo 

ydzGo(x-z)A:e(z-y) = J ds Ks+d^x-- y), 


(59) 


(60) 


which follow from (55) and (59). The parameter e regularizes the singularity in the integration over the proper time 
s at s = 0. We use imaginary frequency, w = in order to avoid inconvenience with powers of i. We use the heat 
kernel regularization in the following way. For the scalar case we take (Ii(£,x) = Kg{x). For the polarizations of the 
electromagnetic case we take (5i(e,x) = PKg{x) which with (10) turns into 


TE : (li(£,x) = 4:Tr^‘^Kg{x.), 

TM : (5i(£,x) = 47r(^2 _ 

P: (5i(£,x) = - d'^)Ks{x.), 


(61) 


and for all cases we take 62 ( 5 , x) = Kg{x.). Now we insert these definitions into eq. (54). For the scalar case we get 

1 

00 = -+/ dse-^^"K,+ 2 ei 0 ) (62) 

9 Jo 
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and for the electromagnetic cases we get 

=-hlTT / ds 

a|| Jo 

1 

4>™ = -hlTT / ds - An)Ks+ 2 e{^) 

“II Jo 

(j>Q = -h 47r / ds - d'^)Ks+ 2 s{^) 

“3 Jo 

Calculated for £ —>• 0, these expressions become divergent according to 

1 1 


x=0’ 


x=0 


(63) 


iscalar _ i 

Z>0 —-r 


+ ^ + 0(v^), 


,TE 

^0 


-TM 

'^0 


/'O 


g 47r 

2 

+ iuY + 


“II 

1 

“II 

1 

0:3 


V 2 Tre 

1 


UJ 


(oyfY/KE^I'^ ‘iyf2nre 

1 


' “ 

2 iijY 


-O(Vi), 


GVStte 3 

where we turned back to real frequencies, ^ = —zw, and give rise to the renormalizations of the couplings, 


(64) 


1 

- M 1 


for scalar 

^ren 

~ g ^ (27r)3/2v^’ 


1 

^ren 

ckll 

_ 1 UJ^ 

0^11 y/2-Ke ’ 


for TE, 

1 

- ^ 1 ^ 


for TM, 

^ren 

0=11 

a\\ 6V27r£3/2 

3V27r£’ 

1 

- 1 1 1 

u? 

for P. 


03 12V2'n-e3/2 

6>/27re ’ 


(65) 


It is seen that all these renormalizations are different one from another and depend in the electromagnetic case on 
frequency. In terms of the renormalized couplings, we get 


iscalar 

^0 


I zw 

jren ’ 


,TE 

*’0 


I 


CHi I 


iTM 

^0 


1 


a 


icu'" 

ren ' ^ 


^0 


1 


2za;'^ 


at 


( 66 ) 


Finally we consider the wave functions which follow from eq. (57) with £ —>■ 0 and (55). In the scalar case we get 


4>(x) = - 




47r((); 


.scalar 


(67) 


repeating the identification (58) and eq. (67) coincides with eq. (14). 
In the electromagnetic cases we get, using (61), 


i^TE(x) = + 


^2 

1 


■^tm(x) = ^ + A||) 

i?p(x) = e*‘^- + -^(w2 + 92)'^ 


x| 

2a;|x| 


( 68 ) 


which are the expressions for the electric fields in the presence of one point dipole with polarizability according to 
(5), provided proper definition of the (/)qS which will be given in the next subsection. 


E. Electrostatic approach 


In electrostatics, also if including retardation, one considers an applied electric field, in our case a plane wave, 

Eappi,(x) =Eo (69) 
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It polarizes the dipole with dipole moment p, located at the origin, whose electric field + V o V) p, together 
with Eappi, (x), adds up to the total electric field, 


E(x) = Eo 


e*kx + (^2 ^ y Q y) 


gi<^|x| 

47r|x| ^ 


The dipole is assumed to be polarized by a local field, Eioc.(x)) according to the second line in (2), 


p = dEioc.(O), 


(70) 


(71) 


at the location of the dipole. Now, in case of a single dipole, the local field is assumed to coincide with the applied 
field at this position, 


Eloc.(O) = Eappi, (0) 


(72) 


and from (70) one comes to 


iui\x\ \ 

j;yi)4E.„,(0), (73) 

This is the electric field including the response of the dipole. It is to be mentioned that in this procedure no 
singularities appear and that the result is perfectly fine and in (73) the polarizability a is the original one from eq. 
(2). For instance, from the second term one can derive the scattering amplitudes of Thomson scattering, Rayleigh 
scattering etc. The properties of the dipole enter through its polarizability and no ambiguity appears. 

Now we compare this approach with the previously considered ones. Inserting for a from (5) and selecting with (8) 
the corresponding modes and components of the electric field, one comes to the conclusion that one needs to include 
the terms with iuj^ in (64) into the renormalization, thus defining 


E(x) = Eappi. (x) + ( (w^ + V o V) 


and one comes to 


1 

- 4 ,ren 

^11 

1 

- 4 ,ren 

^11 

1 

-^ren 


_I , 

0^1 I y/2TT€ 

J_ + _1_ 

1 I 1 

0:3 


3\/27r£ 


12 \/ 27 r 5 '^/^ Q\/2'ks 


ICO 

3 ’ 
2 iuj^ 
3 


for TE, 
for TM, 
for P, 


iTE 

^0 


Oil I 


iTM 

^0 


Oil I 


(74) 


(75) 


in place of (66) and for the 0 ;“'®"- we have to take the original ones appearing in (73) or, with (5) from (2). As 
already mentioned, the renormalization is not unique. In the approach of using an equation like (7) with subsequent 
renormalization of the parameters d, like in the scalar case, these loose their meaning. The renormalized parameters 
a™"-, (65), can be related via equations (68) with the scattering amplitudes like in the scalar case. In the electrostatic 
approach, because of disregarding the self fields, one does not have an equation like (7) but the system (2) with the 
prescription resulting in (73). In this case no renormalization is needed for. It is, if comparing this approach with 
that of subsection IID, that the renormalization (74) is suggested in order to obtain the same result. 

It must be mentioned, that this electrostatic approach is a kind of workaround. It is a known fact, that equations 
(2), considered together, are singular like eq. (7). The workaround is generally accepted and it gives in all known 
cases the correct result. For detail we refer to Chapt. 17 in [12]. 


III. MULTIPLE CENTERS 

In this section we consider a plane lattice of delta functions at the locations 

Xn = an, with n=(ni,n2,0), (76) 

{rii integer) forming a homogeneous lattice in the (x, j/)-plane with spacing a. This lattice has translational invariance 
and, for instance, for the difference between 2 locations 

— ^n—m 


( 77 ) 
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holds. 

We start with formulas with are valid for an arbitrary locations of the centers. For the scalar field we have the 
equation 


A + 5^(5^(x-a„)j <I>(x) =0, (78) 

which was already mentioned in the Introduction, eq. (1). We use the approaches discussed in the preceding section 
for a single center. Eq. (78) is for the scalar field which we consider first. The extension to the electric field is given 
below. 

In the method of selfadjoint extensions we have first to restrict the domain to functions regular the locations of all 
delta functions and then to add the singular modes. In this way, the solutions have expansions around each delta 
function. 


$(x)=crr - - - r + OSE + 0(|x - a„|)V (79) 

V|x-a„| J 

with extension parameter ose- This is equivalent to use boundary conditions (13) at each a„. We chose the same 
parameter asE for all locations not to break the homogeneity of the lattice. 

The solutions for x ^ a„, i.e., outside the centers, read 


<i>(x) = ^ /„ 

n 


gicj|x-a„| 

|x - an| ■ 


(80) 


These may be interpreted a incoming plane wave and outgoing spherical waves from each center. Expanding this 
solution near am, 


$(x) = + /. 


+ y: /n ^ 


^|am—an| 


n^m 




and comparing with (79), we identify 


ikan 


Q^SE = 


“ /m |am - ani 
n^m 


which is a system of equations for the /„, 


iu;|am —an| 


(aSE - *w)/m - V] -r/n = 6 

f ^ a___ — n__ 


zkan 


n^m 




(81) 


(82) 


(83) 


The inhomogeneous solution of this system, being inserted into (80), delivers the solution describing the response 
to the applied plane wave. The homogeneous solution, if existing, describes intrinsic excitations or bound states, 
depending on the application considered. 

Obviously, in the approach with the zero range potentials one comes to the same equations. A similar picture 
appears with the approach with regularization and renormalization which we consider now in more detail in terms of 
separable potentials. For several centers, the kernel (51) takes the form 


Ar(x,y) = (-w^ - A)(5(x -y) + 5 ^ (5i(e, x - an)(52(£, am - y) 

n 


(84) 


and for the Green function one has 

G(x,y) = Go(x-y)-^ f dx'dy'Go(x-x')(5i(e,x'-a„)4>" m'^ 2 (£,am-y')Go(y'-y) (85) 

n.m"^ 


with 


1 

9 


+ 


/ 


n.m 


dxdy 62 ( 8 ,- y)Go(y - x)(Ii(e,x - am). 


( 86 ) 
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Inserting (85) into (52) results in 




which is analog to eq. (83) in the extension approach. 

The generalization of the solution (56) corresponding to an incoming plane wave follows from eq. (85) in the same 
way as before from eq. (56), 

$(x) = e**"’" - ^ f dx' Go{x- x.')Si{e,x' - / dy ^ 2 (e, am - y)e*'"^. (88) 

Taking £ —>■ 0 in the regularized delta functions in this formula, one comes to 

p*<^|x-an| 

<l>(x) = ^ (89) 


Comparison with eq. (80) gives 






establishing the relation to the extension method. 

Starting from here we restrict the discussion to the plane lattice given by eq. (76) and make use of the translational 
invariance (77). In that case, solution (80) has a Bloch wave property. 


As a consequence. 


holds and eq. (83) turns into 


$(x + a„) = e*‘"®"$(x). 


/n = /oe" 


(asE - ^ 


am — a-n | +'ikan 


having simply the solution 


/o = e*' 


asE - iio- di(w, k) 


Here we introduced the notation 


ds(w,k) = 


^ ^za;|an|+ikan 


which is the generic form of the sums appearing in the considered type of problems. For w = 0 and k = 0, Js(0,0) is 
an Epstein zeta function. In (94) the convergence of the sum comes from Qw > 0. 

With (94), (92) and (90), to solution (89) can be written as 


$(x) =e*'^" + /oFb.k(x), 


where we introduced the notation 




exp(fw|x — an| + fkan) 


which is a weighted sum over the spherical waves outgoing from each center. In (96), the coefficient /o can be rewritten 
using (16) and (94), 


1 - / Ji(w,k)’ 
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where / is the scattering amplitude for a single center. The denominator in (98) is clearly the result from multiple 
scattering within the lattice of delta functions. 

Next we consider the corresponding formulas with the regularization in terms of separable potentials following from 
(85) and (86) making now use of the translational invariance (77). For the regularized delta functions we use directly 
the heat kernel (59). In this way we get 


d>(x) = ^ / ds iF,+,(x - 


— 1 — 


(99) 


for the solution (88), where we turned again to imaginary frequencies like in (38). Further we used 


dyKsia^ - y) = 


iky _ — 


( 100 ) 


following from (59). From (86) we get 

in.m = -^ + / dse~ ''^ Ks+2e {an-m), (101) 

9 Jo 

which is, in fact, a function of the difference, in.m = in-m.o- This allows to invert the matrix (j^n.m by Fourier 
transform. Defining 


n 


we get 




-1 g*ka„ 


^ikan 


(j)(k) 


and the solution (98) takes the form 


with 


$(x) = 


47r^(k) 


7^£,a;,k(^) 


7^£.<.,;.k(x) = 47r^ / dse ®''''+**^"iFs+£(x - a„). 


( 102 ) 


(103) 


(104) 


(105) 


This is the sum of the smeared (regularized) outgoing spherical waves from each center. For £ —0 it turns into (97) 
and the relation 


/o 


-1 


(106) 


holds. 

Next we have to consider the renormalization. It is hidden in eq. (101) in the diagonal contributions, 

0n.„=-+/ dse-^«'iF,+2e(0), (107) 

9 Jo 

which is just (j)o in eq. (62). It is obvious that the divergent contributions come only from the diagonal elements 
due to the decrease of the heat kernel in (101) for Un-m ^ 0. In this way, all formulas of section II.3 related to the 
renormalization apply here too. 

Especially for the scalar case we get then in (102) with (65) 


^(k) 


1 

9r 


ICO 

47r 


E 


dse ''^^Ks+2s{an)e 


ikan 


(108) 
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where the prime at the sum means, as usual, to drop the term with n = 0. Here we can put £ = 0 and get 


<^(k) = 


1 

9r 


iuj Ji(a;,k) 

47r 47r 


(109) 


with Ji(a;,k) defined in (95). Inserted into (104) and taken for e = 0, this repeats just (96) with (94) and (46). 

The same can be done for the electric field using the modes corresponding to (10) and the renormalization according 
to eq. (74). Basically, it amounts in inserting the operator P defined in (10), in front of the heat kernel in (108) and 
carrying out the derivatives. Using eq. (61), we get, already for e = 0, 


^™(k) 


0™(k) 


^P(k) 


1 


1 


1 



(47rs)^/^ 

(47rs)^/^ 

(47rs)^/^ 


47r^^ exp(-^ +tkan), 
As 


+ exp(-J+zka„). 


( 110 ) 


Carrying out the integrations over s and using the notations (95) we come to 


0™(k) 

^™(k) 

^P(k) 


1 


CX\ I 


- w^Ji(w,k), 

+ iujJ2{oJ, k) — Jsiuj, k), 

— Ji(w, k) — tw J2(w, k) + Jsiuj, k). 


( 111 ) 


which are the formulas coming for the electric field in place of (109) for the scalar field. 

Now the solutions for the electric field are given by eq. (104) with inserting P from (10) in front of eq. 

(105), which follows from (61) and the substitution of <()(k) by that from (111) for the corresponding cases. Taking 
the limit £ —>■ 0, one comes to the formulas 

£;te(x) = ^ uj^Fu:,k{x) 

(()™(k) 

£^tm(x) = -^‘^>k(x) 

Up(x) = ^ + af) F^,k(x) (112) 

(()^(k) 


where for the </> one needs to insert from (111). We mention that the polarization vector usually appearing in front of 
the plane wave contribution is absorbed in (8) and does not appear here. In fact, Fcj^k(x) with P from (10) applied, 
are the well known formulas for the electric field from a dipole, written here for the specific cases considered. 


IV. SCATTERING ON A TWO DIMENSIONAL LATTICE 


In this section we consider the solutions found in the preceding section for a plane lattice in a scattering setup in 
two cases, for outgoing spherical and plane waves. In both cases, the incoming wave is the same plane wave as before. 
For the spherical setup one considers an outgoing spherical wave at |x| —>■ oo. 


$(x) ~ P'^^ + fs. 

|x|—>-oo 


(113) 


where fsc is the scattering amplitude. For a single center, the solution (14), has already this form since the second 
term in the solution is everywhere, except on the origin, a single outgoing spherical wave and fsc = / holds. 

For multiple centers (76) we have as usual for |x| —>■ oo, 


|x — a„| = |x| — axn + ... 


(114) 
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where x = x/|x| is the scattering direction. From (97) we get 




— mu;xn+iakn 


|x|^oo |x| 


and from (96) and (113) the scattering amplitnde 


fsc = /o X!' 


^ — iaoj'kn-\-iaU.n 


follows. Doing Poisson resummation we get 


fsc = fo 


^Jfwxii -k|| + ^n 


with a two dimensional delta function and the subscript ’||’ denotes the directions parallel to the plane. In this way, 
only scattering in the directions given by 


xii = — kii -n 


is allowed. These are just the v.Laue conditions. Within the zero range potential approach this was pointed out in 
[29]. 

Next we consider an outgoing plane wave. For this it is useful to rewrite FL;,k(a;), eq. (97), using (55), 

/ dn c*p(x—on)+ikan 

11 

Doing again Poisson resummation and carrying out the integrations over pn, one comes to 


7^a;,k(^) — 


The integration over can be carried out too. 


f°° ^ exp(i (k|| + ^n) X|| + ipsz) 
1-00 271"^ — zO 




where we defined 


r„ = yc^2-^k|| + -nj , 

inheriting a positive imaginary part, Sr„ > 0, from uj. Inserted into (96), the solution becomes 

r'„ 


Being interested in scattering solutions not decreasing for \z\ —>■ oo, one needs to have real r„. This implies the 
restriction 

|n| < ^(w-fc||), (124) 

thus restricting the summation in (123) to a finite sum. A restriction to small n appears, for instance, for a small 
lattice spacing a, or, with w = + A:|, for small k^. However, scattering with n = 0 is always possible. For 

the plane wave scattering setup on a homogeneous plane lattice it is useful to introduce Bloch waves in the parallel 
directions. Dividing the wave vector of the incoming plane wave, 

7 27r 

*11 = qH- m, 

" a 


(125) 
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where q is the quasiimpuls, qi < ^ (i = 1,2), we represent the solution in the form 

$(x) = (126) 

n 


and get from (123) 


^n(^) — 


ikzz 


27rt/o 1 

O -r-i ^ 


(127) 


Here m denotes the zone of the incoming wave. The response has scattered waves with, in general, all numbers n. 
From considering this solution for \z\ —>■ oo, it makes sense to define the reflection coefficients 


or, with (106), 


_ 27rz/o 

" a 2 r„ 


—i 

2 a 2 r„^(A:) ’ 


(128) 


(129) 


and Fn (122) is the wave number of the scattered wave moving in z-direction (123). Since we have to account here 
only for functions not decreasing at infinity, the restriction on n to such values for which Fn is real, holds as before. 

Finally we consider the electric field. The structure of the formulas for the scattering into as plane outgoing wave 
remains the same except for the modifications coming from the P’s in (10) and for /o for which we have to use 
/o —>■ —l/(47 r0(k) in (104) with 0(k) from (111). Carrying out the derivatives in P under the sign of the summation, 
we get in place of (123) 


Pte(x) = 
Ptm(x) = 
Pp(x) = 


2 Tri 


E‘ 


girn|z|+i(k|| + ^n)3( 




r„ 


27ri 


q,20TM 

2 m 


^ (k|| + —n)' 


27r + 


a? 4 P{\P) ^ 


^(cc^-F^+ *F„<5(z)) 


r„ 

gir„|z|+i(k||+^n)x|| 

f~r. 


(130) 


The delta function in the last line is the one known for the normal component of the electric field on a double layer. 
It does not enter, of course, the asymptotics for \z\ —>■ ± 00 , which we are interested in. Eq. (130) can simplified a bit 
using (122). 

Introducing in parallel to (126) the corresponding Bloch waves, one comes to the reflection coefficients given by 


r 


TE 

n 


TM 
' n 


r 


p 

n 


—27ri 2 
a2^TE(k)r„‘^ ’ 
-27ri ^ 


— 2 Txi 

a 20 P(k)r„ 


(kll + 



(131) 


which come in place of (128). 


V. THE TRANSITION TO A CONTINUOUS SHEET 

Formally, for vanishing lattice spacing, a —>■ 0, the scattering centers become dense and form a continuous sheet. 
However, this limiting process is quite singular and unexpected results may appear. First of all we mention that in 
the solution (123), for a —>■ 0 only n = 0 delivers a nonvanishing solution. From (122), using = k^i + fc|, we have 
then Fq = fcs and this solution is 


$(x) = ^ ^ ^ik 3 \z\ 


(132) 
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a^fcg ’ 

or r = ro from (128). Here, the factor -^ = p has to be interpreted as density, which must be kept finite. 
Using eq. (98) for /o and, 

‘2t7TZ 

Ji(w,k) = ^-— for a-)> 0, 

which follows from (95) by substituting, in leading order, the summation by integration, we get 


Further we insert from (46), 


q2 _|_ Iq2^2 ’ 


where the second term in the denominator vanishes for a —>■ 0 since in the first term c? jPr must be kept finite. Using 
that in (135), we get 


which is the reflection coefficient tq, (128), resulting from the solution derived in the preceding section for a —>■ 0. 

The transition a —>■ 0 can also be done in the equation (78). Thereby we ignore the problems that the delta function 
in this equation is singular and, thus, that this transition has a rather formal character. The term with the delta 
functions can be written as 


n n 

For a —>■ 0, the summation index n becomes continuous and the sum turns into an integral, which can be carried out, 

- an) - 5 / dn 5(x|| - an) =-^. (139) 

a-s-O J 

n 

As said above, the p= \/a} is interpreted as finite density. With (139), the equation (78) turns into 

(—— A + (/p(5(z))$(x) = 0. (140) 

It has now a one-dimensional delta function potential and is well defined. Taking Fourier transform, 

$k„(z) = f dx||e-*‘^ii"ii$(x), (141) 


the scattering solution is 




These formulas are well known, we use those from the Appendix in [3]. Comparison with (137) shows the reflection 
coefficients are the same provided the relation g = gr holds. This way, taking the limit a —>■ 0 first in the equation, 
or first solving the equation and taking the limit afterwards delivers the same reflection coefficient and it provides 
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another relation for the renormalized coupling constant, and by means of (16), it provides a relation to the parameter 
of the self adjoint extension. 

For the electric field we start from eq. (130) and consider first the TE mode. Again, for a —>■ 0, only the term with 
n = 0 survives at \z\ — >■ oo in the sum and the reflection coefficients (131) become ro —>■ r like in (133). For the TE 
mode we get from the first line in (131) 


with 


which can be rewritten in the form 


„TE 


— 2 'KiLo'^ 
a20o(k)fc3 


00 (k) = 


2 'KiLo'^ 


TE 


1 - 


2'7rza;2ayf"- 


(144) 


(145) 


(146) 


With this reflection coefficient, the scattering solution is the same as (132) with E'^^(x) in place of 4>(x). Equally 
well it can be obtained from (130) for a —>■ 0. 

Like in the scalar case we consider for comparison the transition a —>■ 0 in the equation. 


^-^2 _ ^ _ 4jraiiuj^ - an)^ £^te(x) = 0, 


following from (9) with (10). Doing the same operations as in eqs. (138) and (139), we get 


+ kfi - - 47ra||.;25(z)) 4™(z) = 0, 


(147) 


(148) 


where we also introduced the Fourier transform like in (141). We get the same equation as (140) with the change 
g —>■ —47ra||a;^. This same substitution, as seen, connects (143) with (146). Thus for the TE mode we have, like in the 
scalar case, a well defined transition to the continuous sheet. These formulas coincide with the hydrodynamic model 
and with [3], which is an expected result. 

Now we come to the TM modes for both, parallel and perpendicular polarizabilities. The solutions for a lattice 
with finite a involve 0™(fc) and <j)^{k) which are given by eqs. (111). For a —>■ 0, the lattice sum J 3 {(jj,k), (95), 
converges and behaves as 


Consequently, 


J3(w,k) ~ —. 


a->0 O"’ 


TTM 


(fc) 


a—>^0 Cl 


\k) 


a^O a'-' 


(149) 


(150) 


hold. This behavior is by one power of 1/a too singular. A factor 1/a^ can be absorbed into the density as done before, 
but any additional factor cannot. As a consequence, we do not get any sensible result for a —>■ 0. This singularity was 
pointed out in [2] for the perpendicular polarizability. A somehow unexpected result is that it appears also for one of 
the polarizations (TM) in case of parallel polarizability (which was not considered in [2]). 

Let us consider the transition a —>■ 0 also on the level of the equation (9). For the TM mode, the formal transition 
using eq. (139) cannot be done directly because of the derivatives in A||. We attempt to handle this by Fourier 
transform, 


and get 


dkii 


^TM(x) = y 


+ fcjl + 0z) ^™(^) - 47ra|| “ (*ll + ^n(O) Hz) = 0, 


(151) 


(152) 
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where we used 

/ c^X||e-*‘^ii^ii(c^2_^A||)^,5(x||-an)S™(0)= ^ - (/cn + S™ ^„(z). (153) 

n ^ n ^ ^ 

As expected, the derivatives turned into the corresponding momenta. 

If now doing the transition a —>■ 0, in case the solution decreases for n 7 ^ 0, only the term with n = 0 

survives. This is just what happened for the scalar case and for the TE mode, where it resulted in the well known 
result (140). 

The decrease of ^^(0) can be guessed by solving the equation as in Sect. II.C. In doing so for the TM case, from 

the factor — (A:|| + ^)^) one observes that this decrease is not present for the TM mode. A similar observation 
can be made if doing Fourier transform (151) with the solntion (130). It is again the factor — (A:|| + ^n)^), which 
the TM mode has more, which prevents the decrease. As a consequence, a transition a —>■ 0 in eq. (152) does not 
restrict the summation to n = 0 . 

It is to be mentioned that the restriction to n = 0, if doing it nevertheless, results in the equation 

+ k^i — d\ — 47 rQ;||(u;^ — k^|)( 5 ( 2 :)^ E™(x) = 0 (154) 

for the TM mode. The scattering solution, written in parallel to (132), has a reflection coefficient 

rTM = -^ 2 —, (155) 

27r2Q:| I fes 

where we used = k^| + This equation and this reflection coefficient are there same as in the hydrodynamic 
model for a continuous sheet. 

A similar picture appears for the TM mode in the case of perpendicular polarizability. The equation following 
formally from (9) is 

— A — 47 rQ; 3 (u;^ + 9^) ^ 5(x — a„)^ Ep(x) = 0. (156) 

After doing Fourier transform (151), it becomes 

+ kf^- dfj (z) - iiras = 0 ’ (157) 

n 

Here the problem is in the second derivative of the delta function, 5''{z), which appears in addition to (154). How this 
problem could be handled was shown in [3] by going beyond the dipole approximation. This is, however, not topic of 
the present paper. In addition to 5"{z), eq. (157) has for a —>■ 0 also the problem that the restriction to n = 0 does 
not appear either. Here the missing decrease of A'^||_|_ 22 r^^(0) for k|| —>■ 00 cannot be seen directly from the equation 

as before, however from the solution (130). It has with Tn a nondecreasing factor. 

In this way, in both TM cases, the derivatives, resulting from those in the first line in the right sides of eq. (2), 
prevent a transition to the continuous sheet. 


VI. TRANSITION TO CONTINUOUS SHEET KEEPING A REGULARIZATION 

We have seen that after renormalization, the regularization can be removed, £ —>■ 0, delivering meaningful results. 
These agree with the electrostatic approach where the self fields are excluded. Further we have seen that a part of 
these results becomes singular when the lattice spacing becomes small, a —>■ 0. In this section we ask the question 
what happens if we first take a —>■ 0, keeping £ > 0. The result should be a sheet of finite thickness given by £. Further 
we consider what happens if we take, after a —>■ 0, also £ —>■ 0. 

We start with the scalar case. The solution is given by (104). For Fe_(^_k(x) we use (105) and (59). For a —>■ 0, the 
summation over n turns into an two dimensional integration according to 

^ 5 (on) ^ ^ / dng{n), 


(158) 
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for some function g{n), and we get 
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$(x) = + -=- hJz) e 

2ia^fc3(/)(k) 


zkiix 


ll^ll 


with 


he{z) = 


-iks ds 


Jo "V 4(s + £) 


“TTrTTV - + k||) - 2£kf 


(159) 


(160) 


We accounted for ^ = —zw, which gives 4" ^fl “ ~*^ 3 - The function he{z) describes the ^-dependence on the 
solution to the regularized equation. It is decreasing, 

h,{z) ^ (161) 

| 2|—>-00 


such that the separable regularization gives a sheet of finite thickness without scattering. It is only for £ = 0 that it 
turns into the scattered plane wave, 

ho{z) = . (162) 

This property gives rise to the definition 


1 

2za^fc3^(k) 


(163) 


of a reflection coefficient such that the solution can be written in the form 


$(x) = 


(164) 


In this way, the wave function, if first taking a —>■ 0 and subsequently £ —>■ 0, turns into the same as when doing the 
limits the other way round, eq. (132), with in place of r. 

Actually, the same happens with ^(k). entering (163). From (108) we get with (158) 


^(k)| 
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(165) 


The term present in (108) is subleading for a —>■ 0 and was dropped. If now in (165) taking £ —>• 0, we get 
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which gives in (163) just the same reflection coefficient as given by (137). Thus, for (()(k), the limits commute too. 

Now we consider the TE case for the electric field. Going through the same formulas which resulted in (112), but 
keeping £, we get 

E™(x) = + r™ he(z) (167) 


with 


and 
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TE 
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^™(k) 
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P ds ^-.(C^+kf|)-2£kf| 
a^ Jq ^/s 2s 


(169) 


It can be seen that this integral is finite for £ —>■ 0 and results in the same reflection coefficient (146) as before. Thus, 
also for the TE mode, the limits do commute. 
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For the TM mode we get the same way as before from (112) 

E™{x.) = + r™ h^{z) 


with 
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27rifc3 


(170) 

(171) 


This formula appears in such simple way since the derivatives in A|| act only on 
picture we observe in 0'^'^(k). Here we have from (110), restoring e, for finite a, 

1 1 W /■“ 


I. A slightly more complicated 
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(172) 


As long as e > 0, we can use (158). For £ = 0 we could not do that because of the singularity at s = 0. This is, 
basically, the moment where the limits do not commute. 

Using (158) in (172), we get 


<^™(k)|,^0 = 




ds 


/q + 2e 


„-^i(C"+kf|)-2£kf| 


which has also a finite limit for £ —>■ 0, 


(0™(k)|a=o) 


1 2TTik3 


le^O Q;|| 


(173) 


(174) 


(we remind that the factor is the inverse density and should be considered as finite). Comparison with (150) shows 
now the difference resulting from the orders of the limits. 

The reflection coefficient following from (148) is just the same as (155) and, thus, in the hydrodynamic model, 
provided one makes in (174) the choice a)™' = Oi\\- 

Finally, we consider the TM mode for the perpendicular polarizability. We get from (112) 

F;P(x) = + dl) he{z) £**^11X11 _ (175) 

^ ^ a2A:30P(k) y ^ ^ ^ 

Again, for £ > 0 we have a decreasing function of 2 :. For £ —>■ 0 it turns into 

lim y + heiz) = (k^ - 2(j(z)) (176) 


This way, the delta function in the normal component of the electric field returns. However, it does not influence 
scattering and we can define a scattering coefficient. 
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(177) 


like in the previous cases. The function ^^(k) entering here follows from (110) using (158) 

1 -x/dTr ds Z' ^ 1 
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a™- a2 i/s + 2£ V 2(s + 2£)^ 
However, this function, in opposite to (173), does not have a finite limit for £ —>■ 0, 
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(179) 


In this way, for perpendicular polarizability, taking the limit a —>■ 0 first, does not help. 

It is interesting to remark that, when removing the singular term in (179) by hand, or by including it into a™''', 
insertion of (179) into (177) gives 


rP = 


-1 


1 _ g^fca 


(180) 


which is just the reflection coefficient obtained in [3], eq. (48) (up to a redefinition of a by 47r). This is, given the 
singularities, a rather formal coincidence. 
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VII. CONCLUSIONS 


In the preceding sections we considered all known approaches to a single three dimensional delta function potential 
and a two dimensional lattice of such potentials. There are basically two approaches, the self adjoint extension method 
and regularization/renormalization. Although different, both approaches are related and deliver the same results. We 
reviewed these relations in detail as well as their relation to the ’zero range potential’ method in quantum mechanics. 
As a nice technical tool we reviewed also the regularization in terms of separable potentials. 

The initial motivation for this paper comes from two dimensional polarizable sheets in electrodynamics. There the 
situation is a bit different. Usually one starts from the system (2) of coupled equations for the electric field in the 
presence of polarizable point dipoles and for the dipole moment induced by the electric field in dipole approximation. 
Inserting the dipole moment into the equation for the electric field, an equation with three dimensional delta function 
potential results. In the generally accepted electrostatic approach, as discussed in Sect. HE, the last step is accom¬ 
panied by dropping the self fields of the dipoles. It is shown, this is equivalent to a specific choice of renormalization. 
In fact, by this procedure, the equations alone, i.e., without dropping the self fields, do not form a closed set. As 
discussed in [12], Chapt. 17, this procedure gives the correct result in almost all cases. To our opinion, it is an open 
question how this procedure can be justified as a limiting case from full quantum electrodynamics. 

Using any of the above approaches, the problem of the interaction of the electromagnetic field with a lattice of 
dipoles can be formulated correctly. For example, within a scattering setup, the wave functions and the reflection 
coefficients can be written down (Sect. IV) in the form of known lattice sums, eq. (95). Formulas like (111) and 
(130), (131) allow to calculate any quantity in question, plasmons or the Casimir force between two parallel lattices 
for instance. Also it should be possible to calculate the quantum mechanical bound levels on a lattice generalizing 
[30], where this was done for a one dimensional lattice of zero binding energy. 

In the limit of vanishing lattice spacing, a —>■ 0, one comes formally to a continuous polarizable sheet. However, 
when using a regularization, there is a second limit. For example, in the regularized equation (48) it is e —>■ 0 
and a finite e corresponds to a smearing out of the delta functions. As it is easy to imagine, these limits do not 
commute in general. It is only for a scalar field and for the electromagnetic TE polarization in case of a sheet with 
parallel polarizability, that these limits do commute. For TE polarization, a result comes out which is known from 
the hydrodynamic model for a continuous sheet. For the TM polarization in both cases of polarizability, parallel and 
perpendicular to the sheet, these limits do not commute. 

In the case of first doing the limit e —>■ 0 of removing the regularization first, for the TM polarization the limit 
a —>■ 0 is singular as seen in Eqs. (149) and (150) for the refl;ection coefficients. For the perpendicular polarizability 
this was observed in [2], for the parallel polarizability this is new and it was not expected. One can discuss that the 
reason for this singularity does not depend on the direction of the polarizability but results from the derivatives in 
the couplings in Eq. (10). 

It must be mentioned that for parallel polarizability, the equations turn formally into that of the hydrodynamic 
model. Actually, this is the reason to expect for the TM polarizability a finite result. We showed at which place it is 
impossible to do a —>■ 0 in the equations and where the difference between TE and TM polarizations shows up. 

In the case of first doing the limit a —>■ 0 of vanishing lattice spacing first, for the TM polarization in case of 
parallel polarizability, one comes to a well defined reflection coefficient and the subsequent limit £ —>■ 0 gives the result 
expected from the hydrodynamic model. This result is a bit counterintuitive since here the smeared delta functions 
largely overlap forming a sheet of finite thickness. This sheet is quite specific due to the separable regularization. For 
a sheet of finite thickness with homogeneous permittivity it is known (see, for instance Sect. V in [3]) that the limit of 
thickness to zero must be accompanied by an increase of the polarizability (or permittivity) in order to get a nonzero 
result. In the separable regularization, considered here, this is not the case. 

A similar picture appears for the perpendicular polarizability (here we have only the TM polarization). However, 
in addition here a singularity appears which results from the delta function the normal component of the electric field 
has on a double layer. By removing this singularity by hand, the same result appears which was obtained in [3] by 
stepping back from the dipole approximation. 

In summary, for the electromagnetic field, the results depend on the order of taking the different limits involved. 
We consider this as a hint to go beyond the dipole approximation. 
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